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Matrix product states and their continuous analogues are variational classes of states that capture quantum 
many-body systems or quantum fields with low entanglement; they are at the basis of the density-matrix renor- 
malization group method and continuous variants thereof. In this work we show that, generically, Appoint 
functions of arbitrary operators in discrete and continuous translationally invariant matrix product states are 
completely characterized by the corresponding two- and three-point functions. Aside from having important 
consequences for the structure of correlations in quantum states with low entanglement, this result provides a 
new way of reconstructing unknown states from correlation measurements e.g. for one-dimensional continuous 
systems of cold atoms. We argue that such a relation of correlation functions may help in devising perturbative 
approaches to interacting theories. 



I. INTRODUCTION 

o- 

£NJ Quantum states of many-body systems or fields are character- 
, — ized by their A-point correlation functions. Unsurprisingly, given 
3ieir central status in the respective theories, there are many ways 
which such correlation functions can be book-kept in terms of 
r^as simple as possible mathematical objects. For instance, promi- 
perturbative methods for the description of interacting field 

. theories make extensive use of the relation between high order 

crjr.H two-point correlators fl-4]. These methods, supported by 
Qteerlis' or Wick's Theorem [3,13] (depending on context), give 
+J»ise to a practical way of identifying the propagators as the basic 
Objects for the description of the situation at hand, as well as an 
Siterpretation in terms of virtual processes. 
Q-i In this work we show that, remarkably, generic translationally 
i-Jttvariant matrix product states o-Qj. and their continuous ana- 
logues, cMPS or holographic states Ir9l— Tlltl are completely char- 
r H5cterized by their two- and three-point functions. These states 
j^omprise a variational state class that approximates states with 
^rhmited spatial entanglement well — a ubiquitous property for good 
l/Tjeasons |7|,|g] — and are at the basis of the seminal density-matrix 
^Qfenormalization grou p (DMRG) method |8] and continuous ver- 
j^sions thereof Ir9l 1 1 Qtl - What is more, in our approach states and 
(^corresponding operators can be constructed such that their N- 
£N}oint correlation functions are completely characterized by their 
1 — Correlators of up to arbitrary odd order. We do so by proposing 
^jn explicit construction procedure of how to reconstruct higher- 
• r^rder correlation functions from lower order ones. This insight 
has a number of interesting consequences. 
^ .To start with, a fruitful research program has emerged in recent 
years of revisiting questions in many-body theory within the vari- 
ational set of matrix product states, now seen as a "theoretical lab- 
oratory". This approach has the very appealing feature that some 
links and statements that are in all generality too hard to capture 
analytically can be formulated in a completely rigorous fashion. 
In this mindset, complete classifications of quantum phases have 
been given |[l2tl . new instances of Lieb-Schultz-Mattis theorems 
proven flUl . or phase transitions of arbitrary order identified [14]. 
Our statement provides a new tool to grasp the structure of matrix 
product states and their analogues for quantum fields. 

Our result also identifies matrix product states as a variational 
class that is similar to, but yet beyond, quasi-free approaches. This 
observation may be even more interesting in the light of the fact 
that it is not straightforward to construct natural classes general- 
izing Gaussian states: For example, it is known that any unitary 
evolution generated by quadratic polynomials in the canonical co- 
ordinates maps Gaussian states to Gaussian states. If one looks at 



the closure of the unitaries generated by the quadratic polynomi- 
als and a single further term, say, of third order, one does not 
arrive at a meaningful new variational class, but in fact generates 
all unitaries lfl5ll . 

More practically speaking, our result clearly opens up novel 
ways to think of reconstruction methods for quantum states. Most 
prominently in this direction, is seems unclear how to obtain 
precise knowledge about the quantum states of continuous sys- 
tems such as cold atoms on top of atom chips iflrjll from correla- 
tion measurements — let alone complete tomographic knowledge. 
Such information would, however, be most valuable to analyze 
the experiments made with such systems in and out of equilib- 
rium. The present result shows how continuous-matrix product 
states can be reconstructed from low order correlation functions 
only e.g. obtained from atom detectors lfl7ll . with a phase refer- 
ence obtained by interference. 

Finally, the most important implications may come from the 
description of the fundamental insight into the structure of corre- 
lations as such. In a way, our result shows that states of limited 
entanglement have interesting inherent structural properties, com- 
parable to the structural properties of meromorphic functions. 



II. SETTING 

The main theorem of this paper applies to generic translation- 
ally invariant (continuous) matrix product states in the thermody- 
namic limit. Let us define what this means and fix some basic 
notation. 



A. Matrix product states 

A discrete matrix product state vector of an A-partite spin sys- 
tem with periodic boundary conditions is given by 

|^mps)= Tr[AW[s JV ]...A (1) [s 1 ]l \s N ,..., Sl ), (1) 



where the j4.W[s;] £ <C dxd are finite dimensional complex ma- 
trices. In this work we will focus on the thermodynamic limit, 
i.e. N — >• 00, and the translationally invariant case, i.e. 

[s] for all i, j. The finite bond dimension d will be arbitrary 
but fixed. In this setting, correlation functions of a set of operators 
{Oj } labeled by an index j and with support on (different) sites 
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ik with = ii<...<iAr take the form 



Iq^n) q(in-1 



Tr 



JJV-i 
j\^[jiv]£«iv-ijv 



,.M^E° 



Cf (n), 



(2) 



With MW = Emn^*H ® H 

Oj |n), and the transfer 
matrix E = J2 S ^* H ® ^H- The star indicates complex conju- 
gation of the matrix elements. We have written the distances in a 
compact form as n = (12 — i\ — 1, . . . ,ijv — iN-i — 1) G S^ -1 
and summarized likewise j = (ii, • • • ,iiv). It is possible to con- 
sider finite dimensional and infinite dimensional local systems; in 
the latter case the matrices A[s] have to be chosen such that the 
infinite sums converge. 

The expectation values are invariant under simultaneous con- 
jugation of all M^l and E with some invertible matrix, making 
it possible to consider an equivalent formulation where E is in 
its Jordan normal form (JNF), i.e. E h-> J{E). We call the MPS 
generic if J(E) has non-degenerate diagonal entries fix, . . . , [i^i 
and, moreover, if the largest absolute value occurs only once. We 
order the diagonal elements by their absolute value, in descending 
order. Note that in the thermodynamic limit normalization implies 
Mil 5= 1> where the one with the largest magnitude equals unity, 
i.e. Hi = 1. In the future, we will simply say eigenvectors when 
we mean the right eigenvectors, i.e. E \i) = /Uj \i). The number 
of Schmidt coefficients, and hence the entanglement belonging to 
any contiguous bipartition of regions is limited by 2d. 



B. Continuous MPS 

A one dimensional non-relativistic bosonic quantum field can 
be described in terms of field operators ^>(x) and ^^(x), with 
[*(a;),*(a;') t ] = 5(x - x') and *(x)|0) = 0, where |0) is the 
vacuum. A particular class of one dimensional quantum fields is 
that of continuous MPS (cMPS) or holographic states HOU, 

IV-cMps) = Tr aux [Ve-fo d*Q(*)m+R(*)®*H^ |q) ; (3) 

where Q(x) and R(x) are x-dependent finite-dimensional com- 
plex matrices acting in a d-dimensional auxiliary space. Similar 
to the case of MPS, we focus on translationally invariant cMPS, 
having constant Q and R, in the thermodynamic limit L —> 00. 
It is useful to introduce the Liouvillian (see appendix) matrix de- 
fined as 



T = Q* 



1®Q + R* ®R. 



(4) 



A state of such a quantum field is completely characterized by 
all the possible normal ordered correlation functions of the oper- 
ators ty(x) and W(x) e.g. 



(^(x 2 )^(x 5 )...^(x i )^(x 3 )^(0)). 



(5) 



where the order of position labels is such that they increase in 
size from left to right within the &\ decrease within the "J, and 
= Xi < . . . < Xfq. Correlation functions of cMPS are given 
by expressions involving only the auxiliary space. Let e T °° be 



a short notation for limj,_ 
function, we have 



JL 



In the example of a 2-point 



according to the well-known expressions for expectation values 
in cMPS. Sometimes we might have access to the density only, 
i.e. to the operator ^(x)^f(x), whose correlations are |0 llj, for 
instance, 

(* t (0)* t (x)*(x)*(0)) = Tr [{R* ® R)e Tx (R* ® R)e Toc ] . 

(7) 

Let us go back to Eq. (0. For translationally invariant cMPS, 
we consider the differences between points, Tj = Xi+\ — Xi and 
summarize them in a vector notation r = (ti,T2, • • • , tat_i) G 
R^ 1 . Let the matrices be equal to R*® 1, 1 ® R or R* ® R 
etc. With this notation we represent all A^— th order correlation 
functions in a compact and straightforward way. For example 



Tr 



M [l] e Tr 2M [ 3 ] e 7V 1M [2] e 



Tod 



(3) 



(r) (8) 



(^ t (x)*(0)} = Tr [(1 ® R)e Tx {R* ® l)e 7 



(6) 



with t = (x 2 ,x 3 - x 2 ), j = (1,3,2), and A/W = R* ® 1, 
Ml 2 ! = 1 ® R and = R* ® R. Note that also in this case a 
gauge transformation is possible, corresponding to a simultaneous 
conjugation of T and the matrices by an invertible matrix, so 
that we can always go to a picture where T is in its JNF. The rela- 
tionship between cMPS and channels directly implies 1 1 811 that the 
diagonal elements Ai, A2, . . . , \d 2 of J(T) are closed under con- 
jugation. We call the cMPS generic if J(T) has a non-degenerate 
diagonal and, moreover, the largest real part occurs only once. 
We order the eigenvalues in descending order by their real parts; 
the normalization of cMPS implies that it is non-positive and the 
largest one is zero, i.e. even Ai = 0. Similar to MPS, cMPS have 
limited bipartite entanglement between contiguous regions: the 
number of Schmidt coefficients is limited by 2d. As in the case 
of MPS, we will simply say eigenvectors when we mean the right 
eigenvectors, i.e. e T \i) = e Xi \i). 



III. MAIN RESULT 

In general, to characterize the full state of a quantum system 
one needs to specify all the correlation functions. One may ask 
the following question: "Is it possible to completely characterize 
a (continuous) matrix product state from low order correlation 
functions?" With the only initial assumption of a given bond di- 
mension d, in this work we will show how to: 

1. Certify that the given (continuous) MPS is generic. 

2. Reconstruct the full state of a (continuous) MPS from low 
order correlation functions once the previous property has 
been verified. 

In order to accomplish the first task we need to test if the diag- 
onals of J(E), J(T) are non-degenerate (see the corresponding 
sections). The diagonal elements will appear as poles of the Z- 
respectively Laplace-transform of the correlation functions. The 
second task is more involved since it requires a careful analysis of 
the structure of correlation functions, instead of the mere knowl- 
edge of the diagonal. Both aspects will be studied in detail in the 
following. 



A. Data structure and transformations 

We will use both the Z- and the Laplace transform of correlation 
functions in their multi-dimensional form. For discrete MPS, the 
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Z-transform is applicable 

■2f°(s)= J2 »" 1 ---*S pr Cj W (n) ) Sl! ..., Sjv eC. 

(9) 

Similarly, we have a Laplace transformation of the cMPS correla- 
tion functions 

£f°( s )= [°° d N - 1 Te- s - T d N \r), s u ...,s N e€. (10) 
Jo 

Depending on the correlation data, these transformations will 
not converge everywhere. The key observation is the following. 
A close inspection of the correlation functions reveals that under 
the assumption of a non-degenerate diagonal of J(E), we have 



d- 



fcl,...,fcjV — 1 = 1 



x(^-i) nw - 1 ---(wJ ni ) (ID 



where 



The fact that the diagonal of J(E) and J(T) is non-degenerate 
and has only finitely many elements played a crucial role in the 
derivation of the form of Eq. (TTTT) and (TPfl i. This procedure is 
stable under the influence of noise in the sense that the residues of 
the poles and their positions are continuous in the matrix entries. 
For the purposes of the present work, we will postpone questions 
of the scaling of the computational effort of the classical post- 
processing in the precision of the input as well as in the bond 
dimension of the MPS and cMPS IU9I1 and focus on the possibility 
of reconstruction of higher correlation functions from lower ones 
as such. 



B. Reconstruction theorem 

The form of the Eqs. ( TT3l and dT3b implies that all poles are 
elements of {Ai} and {p" 1 }, respectively. Depending on whether 

the corresponding residues Cj (fci, . . . , fcjv-i) are zero or n °t, 
the transforms of the correlation data may or may not reveal poles 
at these points. This makes it useful to give the following defini- 
tion 



(fcl, . . . , fcjV-l) = 



<1| M^-J |fctf_iXfcjv-i| mV*-^ . . . IhXfcilM^] |1) , (12) 

and {|fc)} denotes the basis where E obtains JNF. Note that in this 
basis E°° = |1)(1|. Considering that |^| < 1 with fix = 1, we 
deduce that the Z-transform — as a function in the complex vari- 
ables {sx, . . . , sjv-i} — converges within the product of unit disks 
around the origin, yielding a meromorphic function. It is possi- 
ble, starting from this region, to reconstruct the whole function 
(its poles and the residues) by analytic continuation. Another way 
of dealing with e.g. experimental data, would be to fit functions of 
the given form using only the unit disk as support. For partial re- 
construction of experimental data, see also Sec. IIV Dl where it is 
shown that poles corresponding to long-range behavior are close 
to the unit circle. Summarizing, we have access to 



fci,...,fcjv-l 



(1 - /i fel si) ■••(!- ^fe„_ 1 sjv-i) 



(13) 



Similar considerations apply to cMPS correlations, which, un- 
der the assumption of non-degenerate J(T), can be expressed as 



fcl, = 1 



x e 



Afc i Tl . . . e Afcjv - 1 ' 



(14) 



with the same symbols cj (fei, . . . , fcjv-i) as defined in Eq. (fT2l . 
The integral in Eq. ( TTOb yields meromorphic functions in higher 
dimensions of the form 



4 w) (s)= E 



(Afc, - si) • • • (A fcjv _ 1 - sjy-i) 



(15) 



According to < 0, the region of convergence of the integral 
in Eq. (TTOb is the product of complex half -planes with positive real 
part. As a meromorphic function, it can be reconstructed using 
this data. Note that now e T °° = |1}(1| and poles reflecting long- 
range behavior now lie close to the imaginary axis. 



Definition 1 (p-number) Given a (continuous) MPS with bond 
dimension d, we define the p-number as the minimum order p 
such that d 2 distinct poles appear in the Z- or Laplace transforms, 
respectively, in at least one correlation function of order less or 
equal to p. If the minimum does not exist we say that the p-number 
is infinite. 

Note that in this definition, we only need one correlation func- 
tion of any subset of operators of interest to show all poles, in 
order to derive the p-number. This provides a solution to the first 
task: If the p-number of a (continuous) MPS is finite, we can di- 
rectly claim that E or T have a non-degenerate Jordan diagonal. 

Now we are going to study in more detail not only the poles but 
the full structure of correlation functions of (continuous) MPS. 
We work in the basis where E or T are in JNF. If a function Z> 



or £j is given, each coefficient cj"' 1 (k±, . . . , fcjv-i) can be ex- 
tracted by finding the residue of the corresponding multi-pole. We 
can finally state the main theorem: 

Theorem 1 (Computing higher from lower correiation functions) 

A generic, translationally invariant (continuous) MPS in the 
thermodynamic limit with p-number p is completely characterized 
by the correlation functions of order I < 2p — 1. 

Proof: We merely need to consider the case in which p is finite. 
From the definition of the p-number we know that J(E), respec- 
tively J(T), have a non-degenerate diagonal in the JNF, whose 
entries can be recovered from (the pole structure of the transforms 
of) the correlation functions of order £' < p. This reduces the re- 
construction of correlation functions to the reconstruction of the 
coefficients in Eq. (TT2b . The problem is now to express every co- 
efficient 



(fcl 



,k N - 1 ) = M^: 1 \M^ N - 2] 



(16) 



each associated with a unique set of poles, in terms of low order 
coefficients with £ < 2p — 1. From the definition of the 
p-number we know that for each index k there is at least one non- 
zero coefficient cj.?j\ (..., fc, ...) ^ 0, with p(k) < p and fixed 
j(fc), having k as one of its indices. This allows us to write d 2 
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different versions of the identity (this is a scalar) 



71/ 



>(*)) 
c j(fe) 



(...,*,...) 



(17) 

where fc = 1, 2, d 2 , and the symbol * stands for other indices 
which are irrelevant. Now we reorder the matrix elements in the 
numerator by shifting all matrix elements on the r.h.s. of the index 
k simultaneously to the l.h.s., leaving the order of the other indices 
untouched, i.e. in the following way 



M \m t 



. M 



>(*))/■ k ^ 

(18) 

We can finally put all these resolutions of the identity between the 
matrix elements in Eq. ( TTol l 



(&!,. . . , fcjv-l) 



M [i«-i] 1 (fc JV _ 1 ) M b«-2] 1 (fe JV _ 2 ) jWwH (19) 

l,fe w _i fcjv_i,fejv— 2 Kill v y 

We recognize, in the numerator, several new strings of matrix el- 
ements resulting from the insertion. They have the same structure 
as in Eq. (fT~6T > but a lower order I < 2p — 1. This means that, for 
every N, all the coefficients cS N ^ can be written in terms of 
with I < 2p — 1. In other words, correlation functions of order 
less or equal to 2p — 1 are enough to reconstruct all the others. 
This proves the validity of the theorem. □ 



IV. DISCUSSION 
A. Example 

It is instructive to consider the following case. Given an MPS or 
a cMPS with finite d, let the operators Oj and the state be such that 
the corresponding matrices = XM^X~ X have only non- 
zero elements. Here J(-) = X ■ X^ 1 is the conjugation that takes 
E, T to their JNF. Note that the probability to have this situation 
in an experiment, or using a randomized (continuous) MPS and 
operators Oj, is one. Under this condition, all two-point function 
transforms show all the poles, hence p = 1. Computationally, all 
residues of all the poles of all iV-point functions with N < 3 can 
be obtained. Hence we can, using the construction above, give 
explicit formulas that express all TV -point functions in terms of 
the 2- and 3-point functions. See also Section ITV B I 



B. Applications in tomography 

The framework established here opens up ways to reconstruct 
unknown states from correlation data alone, in the sense of tomo- 
graphic reconstruction. Such tomographic reconstructions based 
on data sets such as MPS and cMPS are important not only be- 
cause MPS (and cMPS) form such a ubiquitous and important 
class of quantum states, as the states that capture situations of 
low entanglement well. For scalable reconstructions of states of 
quantum many-body systems, one necessarily has to make use 
of data structures that are efficient in the system size. In this 
way, the approach proposed here gives rise to a complementary 
picture of viewing tomography of MPS to the method of Ref. 



112 1H . where the reconstruction is based on tomographic estima- 
tion of certain reduced states. The ideas presented here seem more 
promising, however, for the tomographic reconstruction of quan- 
tum field states and continuous systems, for which no method is 
known altogether. Notably, from correlation data of atoms count- 
ing experiments Il6lll7ll with or without a phase reference, a con- 
struction of quantum field states can be done along the above lines 
JH. 

The aim of a reconstruction is to give (a representation of) a 
positive functional on (a representation) of operators correspond- 
ing to our measurements that generates the given data. This func- 
tional itself will be given by formulas [2] and [8] respectively, hence 
the quantities to be derived from the data are matrices J(E), J(T) 
and M^s. Because the origin of the formulas of our recon- 
struction is a (continuous) MPS ansatz-modified by some gauges 
only-we are already certain of the positivity of such formulas. We 
consider two more aspects. 

First, any reconstruction procedure can at best find a represen- 
tative of the equivalence class of all (continuous) matrix product 
states that are identical up to a gauge transformation. From Eqs. 
© and (O, we immediately derive that simultaneous conjuga- 
tions of the transfer matrix E and the operator matrices — 
analogously the Liouvillian T and M^l — by an invertible ma- 
trix leaves the expectation values invariant. So does the trans- 
formation Q i-> Q + at, which changes the norm of the cMPS. 
These are gauge transformations and do not change the observable 
physics. Some of these transformations have been performed to 
obtain the pole structure in the explicit form, diagonalizing the E 
and T and ordering the diagonal elements. Hence we have "fixed 
the gauge", but only partially. 

Second, we need a promise of the dimensionality d of the ma- 
trices. The p-number then tells us if the experimental data is con- 
sistent with this promise and if we "see enough structure in the 
data", i.e. if we can use it for a full reconstruction of one represen- 
tative. As we have seen in Section ( 1IV At in all cases except a set 
of measure zero, the matrices contain no zeros as elements 
and the p-number is one. This allows for a reconstruction along 
the following lines. Let us assume for simplicity to have only one 
operator M. The case of more is then straightforward. 

After the observed data has been processed, we have obtained 
the values of the variables occuring in Eqs. dTTb (or Eq. (TBI)). 
We only need the two- and three-point functions. In particular, 
we found the position of the poles {/i^ 1 } (or {Xi}) of the trans- 
forms of the correlators. The diagonal elements of the matrix 
J(E) (or J(T)) are in one-to-one correspondence to the poles 
of the two-point function when p = 1, We also found the residues 
cj (fci , . . . , fcjv-i), for N = 2, 3, which determine the operator 
matrices M. We now need to reconstruct M and we are done. 

If one uses a scalar ^ 



(20) 



then the Ms change, but the two- and three-point functions (and 
the pole structure) do not. This gauge comes from the conjugation 
of E (or T) and the M s with an invertible diagonal matrix and 
has not been fixed before. We use this freedom to set M\ t k = 1. 
Accordingly, we obtain, via Eq. ( fT9] l, directly from the residues of 
the poles 

M kA =cW(k). (21) 
The remaining elements of M are obtained from Eq. (fT~9b via 

M M = C ( 3 )(M)/c (2) M. (22) 
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This is enough for the reconstruction. If we have more data (e.g. 
higher iV-point functions), it can be used as a consistency crite- 
rion for the assumptions. If the data is noisy, the formulas will not 
be consistent, but maximum-likelihood estimations can be per- 
formed and can profit from more data and the relations discussed 
in this article 12211 . 

C. Arbitrary odd order dependence 

To construct (continuous) MPS together with a set of operators 
that is dependent on correlators up to an arbitrary odd order, in- 
stead of two- and three-point functions, one has to solve for the 
equations given in Eq. ( fT2b . The aim is to choose J(E), respec- 
tively J(T), and the operators {Oj} such that they do not fulfill 
the p— number condition up to a desired order p. Essentially, the 
corresponding matrices {M^} have to be sparse enough. As the 
matrix dimension d is to be chosen freely, this requirement can be 
met in a construction. 



D. Judging the importance of poles and the treatment of infinite 
dimensional auxiliary space 

The decay behavior of contributions to the correlations follows 
directly a) in the case of MPS from the absolute value of the cor- 
responding poles, and b) in the case of the cMPS from their real 
part. (Note the relation between poles and diagonal elements: 
Si = /i" 1 and Si = Ai, respectively.) Hence their position on 
the complex plane is decisive. This means that the MPS poles de- 
scribing slow decay are situated close to the unit circle, and the 
cMPS poles of this kind are close to the imaginary axis. 

As interesting states of most physical models in ID have an 
infinite Schmidt number even if they obey an area law and have 
quickly decaying Schmidt coefficients Q, the order of the poles 
corresponding to their decay behavior is important when recon- 
structing only a subspace of the auxiliary space. According to 
the consideration above, it is hence feasible to compute only rel- 
evant MPS matrix dimensions from a given set of experimental 
correlation data, retaining e.g. only long-range behavior in the de- 
scription. 

V. OUTLOOK 

A stimulating insight is given by the mathematical structure 
of the correlations, which are, as shown, related to meromorphic 
functions with interdependent residues. The structure of correla- 
tions is moreover linked to the quantitative limitation of entangle- 
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ment between spatial regions on a fundamental level. A possible 
deeper significance of this formal assumption remains to be found, 
but similar limitations, such as the Bekenstein bound, come to 
mind. 

One might speculate that based on the findings above — in par- 
ticular the relations of higher order correlations to two- and three- 
point functions — some new insight into diagrammatic perturba- 
tive methods could be obtained. One future direction of inves- 
tigation are virtual processes. As mentioned in the introduction, 
in well-known perturbative methods for interacting field theories, 
the states of the interacting theory are described in terms of states 
of the non-interacting counterpart. This leads to a description in 
terms of quasi-free states, which are determined by their two-point 
correlators, and puts the focus on propagators and an interpreta- 
tion of the theory in terms of virtual processes. While this ap- 
proach allows to predict experimental data with very high pre- 
cision [21], it is not well understood on a fundamental level, and 
gives rise to conceptual problems l24ll . Hence, transfer of the 
structures discussed in this paper to a relativistic setting might be 
interesting. We also note a similarity of the transforms of (con- 
tinuous) MPS correlations with expressions from conformal field 
theory. 

Of course, the relations underlying our approach are already 
summarized in a computationally efficient formal framework: the 
family of states known under the name of MPS, cMPS, and in- 
stances of PEPS. This means, given experimental data with the 
relations above, an optimal book-keeping device would be, e.g. 
, an MPS. However, say, series expansions in a perturbation ap- 
proach starting from MPS and potentially leaving this class of 
states due to closing gaps etc. give rise to different sets of terms 
with different structure, and the MPS scheme is only one possi- 
ble way to interpret it. A different summation order might yield a 
different optimal book-keeping device in this context. 

In a related line of thinking, there is the possibility that the 
meromorphic structure of the correlations, together with the inter- 
dependencies of the poles, enable a different mathematical under- 
standing of the underlying renormalization procedure. Such an 
understanding might help to find variations of the renormalization 
procedure, including possibly meaningful and computationally ef- 
ficient extensions of MPS and cMPS to higher dimensions. 
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Using the isomorphism relating a density matrix defined on a d- 
dimensional Hilbert space to its vector form in a d x d-dimensional 
Hilbert space, 

d 

p{x) h+ \p(x)) = PM\i)\i), (25) 
the above generator takes the form 



L = -iH* ® 1 +il ® H 

-\Yj ( R J R j 1 + 1 ® R ] R i ~ 2R *j ® R i) ■ < 26) 
3 



The Liouvillians encountered above are Liouvillians in this sense, 
where we merely have a single Lindblad operator R. With 



VII. APPENDIX 

A. Quantum dynamical semi-groups 

There is a strong formal relationship between cMPS with bond 
dimension d and quantum channels which are elements of a 
completely positive continuous one-parameter semi-group in d- 
dimensional Hilbert spaces, called Markovian quantum channels. 
For such elements T x there exists a generator C with £*(1) = 
such that 

T x = e xC (23) 

for all x > 0. The generator L of the family of quantum channels 
is necessarily of the form 

C(j>) = -i[p, { R ) R iP + P R ) R i - 2R 3P R ]) ■ ( 24 ) 





Q = iH - -R^R, (27) 

one has 



L = Q*®1 + 1®Q + R*®R, (28) 
which is Eq. Since the spectrum of £ is taken to be non- 
degenerate, there exists a unique stationary state, given by 

Pss = lim e xC (p) (29) 



for every p. The generator is also gapped, so that convergence 
to the stationary state is exponential. Such families of quantum 
channels are sometimes referred to as being relaxing. 



